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ABSTRACT 


This  report  establishes  rigorous  stability  criteria  in  terms  of 
transfer  matrices  and  in  addition  undertakes  a  careful  and  critical  ex¬ 
amination  of  the  concepts  of  complete  controllability  and  observability 
as  applied  to  the  closed  loop.  In  particular  it  is  empahsized  that  the 
latter  properties  are  neither  necessary  nor  desirable  attributes  of  a 
practical  multivariable  feedback  system. 


MODERN  CLASSICAL  MULTIVARIABLE  FEEDBACK  CONTROL  THEORY-FART  I 


1.  Introduction  and  Notation.  The  objective  of  thi3  report  (which  we  project 
shall  be  issued  in  four  parts)  is  to  rewrite,  extend  and  enlarge  the  relatively  new  field 
of  classical  multivariable  feedback  control  in  the  light  of  the  gains  and  insights  pro¬ 
vided  by  modern  optima)  control  theory.  Part  I  establishes  rigorous  stability  criteria 
in  terms  of  transfer  matrices  and  in  addition  undertakes  a  careful  and  critic?)  ov=»mina>- 
tion  of  the  concepts  of  complete  controllability  and  observability  as  applied  to  the  closed 
loop.  In  particular  it  is  emphasized  that  the  latter  properties  are  neither  necessary 
nor  desirable  attributes  of  a  practical  multivariable  feedback  system.  Part  II  (in 
preparation)  is  devoted  to  a  new  and  general  theory  of  broadband  plant  equalization  via 
feedback.  More  specifically  the  following  problem  is  poBed  and  solved:  Given  the 
transfer  matrix  P(s)  of  a  plant  with  asymptotically  stable  "hidden"  modes  what  are  the 
necessary  and  sufficient  conditions  for  the  matrix  T(s)  to  be  realizable  as  the  input-  • 
output  transfer  description  of  a  dynamical  asymptotically  stable  closed  loop  incorporating 
the  plant  and  an  asymptotically  stable  controller  and  feedback  sensor?  In  general  the 
restrictions  on  T(s)  are  a  direct  consequence  of  the  non-minimum-phase  character  of 
P(s)  and  may  be  expressed  as  generalized  gain- bandwidth  constraints.  These  constraints 
reveal  the  possible  tradeoffs  between  dc  gain  (steady-state  error)  and  transient  response 
(overshoot,  etcetera).  In  Part  III  the  topology  is  delimited  to  conform  with  that  in 
current  use  in  present-day  process  control.  For  example,  it  is  technologically  expedi¬ 
tious  to  assign  to  each  controlled  variable  its  own  controlling  loop.  This  of  course  in 
no  way  implies  a  non-interactiig  control  scheme.  In  fact  this  latter  technique  is  shown 
to  have  some  serious  deficiencies  both  from  a  practical  and  theoretical  point  of  view. 

In  part  IV  we  undertake  to  extend  the  above  results  to  plants  and  feedback  control 
arrangements  with  intrinsic  time  lags.  Particular  emphasis  is  attached  to  the  problem 
of  anolgue  simulation  which  at  the  present  writing  appears  to  be  susceptible  to  a  direct 
and  decisive  solution.  These  considerations  lead  us  naturally  to  the  study  of  the  decom¬ 
position  theory  of  multivariable  rational  matrices. 

As  regards  notation  most  of  it  is  self-explanatory  but  nevertheless  we  point  out 

that  n-dimensional  column-vectors  ore  written  x,a  or  ir.  the  alternative  fashion 

x  =  (x,,  x,,  .  .  .  ,  x  )'  whenever  it  is  desirable  to  exhibit  the  components  explicitly.  More- 

over,  A  and  dot  A  stand  for  the  transpose  of  A  and  the  determinant  of  A,  respectively. 

Rv  O  ,  O  ard  1  we  denote,  ir.  the  same  order,  the  mxn  zero  matrix,  the  n -dimen- 
m, n  ~u  n 

sionnl  zero  c  .lun  n-vector  and  the  nxn  identity  matrix. 

Although  most  of  the  results  in  nart  T  ha  "e  been  known  Co  the  author  foi  sc 
vvr.r  many  of  them  have  not  been  published  by  him  prior  to  tin  appearance  cf  thi  rejv  rt 
i  is  <  ha'  e  been  rat!  .  r  careful  in  the  matter  of  credits  and  if  Ref.  I  appears  next  * 
li  i  1  ,  means  Ui.it  to  the  best  of  the  author1  ■  knowledge  the  theert  n:  •.nd  its  ,  . 
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have  a  pi  nrfed  fry  the  first  bint  in  Ref.  Z. 


2.  The  Basic  Theorems.  A  schematic  of  the  general  feedback  control 
arrangement  (GFCA)  studied  in  this  paper  is  shown  in  Fig.  1. 
comparator 


Fig.  1.  The  General  Feedback  Control  Arrangement  (GFCA) 


Assumption  1  (A^).  The  controller,  plant  and  feedback  sensor  are  assumed  to 
be  linear  time-invariant  dynamical  systems  possessing  the  respective  state-variable 
descriptions 


yoTcScw+oc«e(«) 

(1) 

IcW=Hc2c,t,  +  Jc2c|t)  : 

(2) 

ip(t)  =  Fp5!>(t)  +  Gp2p(tl 

(3) 

y_(t)  =  H  x  (t)  +  J  u  (t)  ; 

~P  P~P  P~P 

(4) 

xf(t)  =F^xf(t)  +  Gfuf(t) 

(5) 

yf  (t)  -  HfXf(t)  +  JfU^t)  . 

(t>) 

As  usual,  u,x,y  are,  in  the  same  order,  generic  designations  for  "input",  "state"  and 
"output".  Of  course,  in  the  situation  depicted  in  Fig.  1,  u(t)  is  the  actual  GFCA 
reference  input,  y(t)  the  corresponding  output  and  e(t)  the  feedback  error.  It  follows 
from  Fqs.  (l)-(6)  that  the  associated  transfer  matrices  are  given  by 

C(s)=J  +H  (si  -F  }"^G  =  m  x  r  matrix;  (7) 

c  c  v  c  c 
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and 


where, 


P(s)=J  +H  (si  -F  )"^G  =  n x m  matrix  (8) 

P  P  vp  P  P 


F(s)  =  J^+H^sl^  -Fj)  ^G^  r  xn  matrix,  (9) 


Vc=  dimensionality  of  controller  =  size  of  Fc; 

V  =  dimensionality  of  plant  =  size  of  F  ; 

P  P 

dimensionality  of  feedback  sensor  =  size  of  F^. 


Clearly,  because  of  the  particular  interconnection  of  the  blocks, 

uc(t)  =e(t)  =  uUJ-y^t)  , 

up(t)=yc(t)  * 

uf(t)  =  y(t)  . 


(10) 

(ID 

(12) 


The  sizes  of  the  column-vector  functions  of  time  u(t),  e(t),  Up(t),  y(t)  and  y^(t)  are 
rxl,  rxl,  mxl,  nxlandrxl,  respectively. 

Let  the  Laplace  transform  of  a(t)  be  denoted  by  a(s).  Then,  by  a  straight- 
*  ~ 
forward  analysis  we  find  that 


and 


where 


y(s)  =  T(s)u(s) 


e(s)  =  E(  s)u(s) 


and 


T(s)  =  [ln+  P(s)C(s)F(s)r1P(s)C(s) 
=  P(s)C(s)[l  +  F(s)P(s)C(s)] 


E(s)  =  [lr+F(s)P(s)C(s)]_1  . 


■1 


(13) 

(13a) 

(14) 

(15) 

(16) 


T(s)  is  n>:<  and  E(s)  is  rxr. 


•wBWWJnrrw 


Observe  that  (14)-(  16}  make  sense  i£  and  only  if 

det[lr+F(s)P(s)C(s)]*  0  .  (17) 

Conversely,  if  condition  (17)  is  satisfied,  then  T(s),  as  defined  by  (15),  is  the  unique 
nxr  rational  matrix  function  of  the  complex  variable  s  ”o  + jU)  relating  y(s)  to  u(s)  in 
the  manner  prescribed  by  (13).  Furthermore,  it  should  be  evident  that  if  (17)  is 
violated  e(t),  yc(t),  y(t)  and  y^(t)  are  not  uniquely  determined  by  the  input  u(t)  for  a 
specified  initial  state  xc(0),  £^(0),  x^{0).  Now  even  granting  (17)  it  does  not  follow  that 
the  GFCA  of  Fig.  1  is  dynamical,  i.  e.  ,  it  is  not  necessarily  true  that 

limit  T(s)  =  finite  matrix  .  (18) 

s  -*  • 


The  possession  of  a  pole  by  T(s)  at  s  =  «  means  that  GFCA  is  performing  differentiation 
as  well  as  integration.  More  to  the  point,  suppose  the  input  is  composed  of  r  Heaviside 
unit  step  functions  u  j(t)  of  arbitrary  amplitudes: 

u(t)  —  u  u  (t)  , 

~  ~o  -  1 

uQ  an  arbitrary  constant  rxl  vector.  Then, 

y(e)=i  •  T(s)uo 

and  invoking  the  Laplace  transform  initial  value  theorem  it  is  seen  that  y{0+)  is  well- 

defined  for  every  choice  of  u  only  if  (18)  is  valid.  If  this  is  indeed  the  case, 

~o 

y(0+)=T(»)uo  . 

Theorem  1.  Under  assumption  A,,  the  GFCA  shown  in  Fig.  1  is  dynamical  if 

if  1 

and  only  if 

0  ^  det[l  +F(»)P(»)C(®)]  =det[l  +  P(»)C(»)F(«)]  .  (19) 

r  n 

Or  equivalently,  in  view  of  the  identifications 

F(°°)  ~  Jj  , 

P(»)  =  Jp,  (20) 

CH=JC. 

_ 

lx  A  is  nxr  and  B  is  rxn,  det(l  +AB)  =  det(l  +BA). 

n  r 
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if  and  only  if 


0  t  det(  1  +  J-J  J  )  =  det(  1  +  J  J  J.)  .  (21) 

rfpc  npcf 

Proof.  Suppose  (19)  holds.  Then  (17)  it  certainly  true  and  T(s),  as  defined  by 
(15),  is  meaningful.  Moreover,  (19)  and  taken  together  imply  that  noth 

Cl^FisJPfsXZU)]-1 

and  T(s)  are  regular  at  s=“.  Inversely,  if  GFCA  is  dynamical  then  of  necessity  y(t) 
is  uniquely  determined  by  u(t)  and  the  initial  state.  Hence  (17)  must  hold  and  T(s),  as 
given  by  (15),  makes  sense.  In  addition,  T(®)=  finite  matrix  whence,  because  of  assump¬ 
tion  1,  F(s)T(s)  is  also  regular  at  infinity.  But 

F(s)T(3)  =  lr-[ir+F(s)P(3)C(s)]'1  (22) 


and  consequently 


det[l  +  F(®)P<®)C(®)]  *  0  ,  Q.  E.  D. 
r 


Motivated  by  theorem  1  we  introduce  assumption  2. 

Assumption  2(A^).  The  GFCA  of  Fig.  1  satisfies  the  constraint 

0  f-  det[l  +  F(«)P(®)C{»)]  =det(l  +  J.J  J  )  .  (23) 

r  ripe 

The  true  significance  of  A^  is  perhaps  best  brought  out  by  an  examination  of 

the  state  equations  of  GFCA  when  exhibited  in  terms  of  the  coefficient  matrices  F 

n  c,p,  i 

G  .,  H  ,,  J  ,  describing  the  dynamical  behaviour  of  controller,  plant  and 
c , p, I  c, p, I  c,p,z 

feedback  sensor.  Put 


=  (v  +v  +v.)xl  vector  . 
c  p  i 


£o(t)  * 


*  W 


Then,  after  some  straightforward  algebra  we  find  that 

xQ(t)  “Fo~o(t)  +Go-(t>  ,  t24' 

y(t)  =Hoxo(t)+Jou(t)  , 


D 


(25) 


ear  in  mind  that 


0  i  det(l  +  J  J  J.) 
n  pc  f 


det(  1  +  J,J  J  )  =  det(  1  +  J  J,J  ) 

rfpc  mcfp 


I 


H  >=(1  +J  Jf>_1[H  l-J  HJ;  '  =(1  +J  J-)-1J 
o  npf  L  p*  p  f*  d  '  n  p  f  p 


L  f  pj 


Although  (26)  looks  very  formidable  it  may  be  recast  in  a  remarkably  simple 
form.  By  inspection, 


F  =  (F  +  F  +F  )-(G  +  G  +  G,)L" 

o  c  p  I  c  p  I 


JfJ  J.  1 
{ P  f  r  . 
TTTTTJ  (H  +H  +  H-)  , 
m  c  f  c  '  c  p  V  ’ 

Tr—rr^trr 

p  n  p  c 


where 


L  =  ( 1  +  J  J  J  )  +  ( 1  +  JJ.J)  +  ( 1  +  J  J  J.)  . 
ripe  mefp  npcf 


With  the  aid  of  the  identity 


J,J  J,  1  .  -I  O 

fpfr  c  in  m,n 

-1  T7TJ  O  J  -1 

mefe  n,  rp  n 

rj— tt— rr  i  o  j, 

p  npc  r  r,  mf 


we  can  rewrite  (35)  as 


F  =(F  +F  +F,)-(G  +G  +G,)  O  J  -1  (H+H+H.). 

o  cpf  cpf  n,  rp  n  'c  p  f 

10 
r  r,  m  i 


To  determine  the  stability  properties  of  GFCA  we  must  find  the  eigenvalues 
of  F  .  That  is,  it  is  necessary  to  fix  the  location  of  the  zeros  of  the  polynomial 


^  Q(s)  h  det(sl  -F  ) ;  i.  e.  ,  the  roots  of 
o 


det(sl  -F)  =  0,  v=v+v+v,  . 
V  o  o  c  p  f 

o  ^ 


According  to  a  well  known  result  in  the  theory  of  partitioned  matrices, 
*  A  -i  B  denotes  the  "direct  sum"  of  the  two  matrices  A,  B. 
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(39? 


i 


dct 


‘11 


*21 


‘12 


22 


=  det(An-AI2A-2A21)dctA22 


=  det  (A22-A21A11AJ2)det 


provided  the  inverses  exists.  Hence, 


Ao(s)=det(slv  -F)=±tf  *.deti 
o 


- 

« 

J 

-1 

O 

H 

O 

O 

c 

m 

m,  n 

c 

m,vp 

m,  v{ 

O 

J 

-1 

O  ,  v 

H 

O 

n,  r 

P 

n 

n  c 

P 

n,vf 

1 

O 

J- 

0  ,  v 

o 

H, 

r 

r,  m 

f 

r’  c 

„r'vp 

£ 

c 

“13 

0 

P  -si 

0 

0 

c 

v  .  m 

y.  n 

C  V, 

v.,  v 

V  ,  V, 

c 

c 

c 

c  D 

c  f 

O 

G 

O 

0 

F  -si 

O 

\L*  r 

p 

v  ,  n 

V  ,  V 

P 

P 

P 

P  c 

r  P 

p  { 

O 

O 

G, 

O 

o 

F,-  si. 

v  ,  r 

v  ,  m 

f 

V.,  V 

V,,  V 

f  \>£ 

P 

P 

I  C 

f  P 

where 


L  TJ  J  J  I 

n  p  c  f 


rj 

-1 

O 

c 

m 

m,  n 

=  det 

O 

J 

-1 

n,  r 

P 

m 

1 

O 

J, 

r 

r,  m 

f  J 

Using  (39)  with  ^22=*f"E*v.  an^  recognizing  that 

F(s)=Jf+Hf(slv-Ff)-1Gf 

we  get  (all  subscripts  on  the  zero  matrices  are  now  dropped) 


(39a) 


(40) 


(40a) 


Aq(s)  =  J:T|  *  '  det 


J  -1  OH 

cm  c 

O  J  -1  O 

p  n 


1  O 
r 


F(s)  O 


O 

H 

F 

O 


GOOF  -si  O 
C  C  V 

c 


O  G  O 
P 


F  -si 

P  v 


Af(s)  , 


(41) 


Aj-(s)  =  det  ( s  1  ^  -Ff) 
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I 

i 

1 

I 


t 


Repeating  the  process  two  more  times  in  succession  with  A^-F  -sl^  and 


A^sF^d^  ,  respectively,  one  obtains 
c 


Ao(s)=±n  *’det 


C(s) 

-1 

m 

O 

O 

1 

r 

P(b) 

0 

-1 

I\ 

F(c) 

_ 

•  A£(b)Ap(s)Ac(b)  , 

(42) 

A  (s)=det(slv  -F);  Ac(s)  =  det(slv  -Fc)  . 
r  p  ^  c 


To  calculate  the  determinant  appearing  in  (42)  we  simply  interchange  block  columns  to 
form  the  matrix 


1 

O  C(s) 

rn 

P(s) 

-]  O 

11 

O 

F(s)  1  | 

(43) 


and  then  apply  (39)  two  last  times  to  reach  the  very  informative  formula 


Ao(s)=n_1-detCln+F(s)C(s)F(8)]-Ac(B)Ap(B)Af(s)  ,  (44) 


ti  =  det(ln+ JpJcJf)  =  det  [ln+P(®)C(®)F(®)] 


The  result  embodied  in  (44)  is  undoubtedly  one  of  the  most  important  in  the  entire  field 
of  classical  multivariable  control.  Our  next  theorem  is  immediate. 

Theorem  2  [2,  3],  Under  assumptions  A^,  A^,  GFCA  is  asymptotically  stable  if 
and  only  if  the  scalar  function 


det[ln+P(s)C(s)F(s)]- Ac(s)Ap(s)A{(s)  (45) 

is  devoid  of  zeros  in  Re  s  >0. 

Unfortunately  the  test  for  stability  provided  by  theorem  2  cannot  be  carried  out 
n  terms  of  transfer  matrices  alone  but  requires  knowledge  of  the  three  polynomials 
Ac(s),  Ap(s)>  A^s),  quantities  which  depend  on  the  internal  structure  of  GFCA.  At  the 

The  sign  ambiguity  is  removed  by  comparing  coefficients  of  the  highest  powers  of  s. 


9 
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present  level  of  generality  nothing  can  be  done  to  remedy  this  situation  .  However,  if 
one  additional  constraint  it  imposed  on  the  controller,  plant  and  feedback  sensor  it  is 
possible  to  eliminate  this  defect  completely. 

Assumption  3  (A.^).  The  controller,  plant  and  feedback  sensor  are  completely 
controllable  and  observable  realizations  of  their  respective  transfer  matrices  C(s), 
P(b).F(s). 

The  following  algebraic  characterization  of  A  is  well  known  [4].  Let  $  (s), 

**  j  c 

ifp(s)  and  f^s)  be  the  characteristic  denominators  of  C(s),  P(s)  and  F{s),  respectively. 
Then,  A^  is  valid  if  and  only  if 

AJb)  =  » 

c  c 

A(a) = i  (s)  ,  (46) 

P  P 

Af(s)  =tf(s)  . 

Theorem  3  [2,  3],  Under  assumptions  A^-A^,  GFCA  is  asymptotically  stable 
if  and  only  if  the  scalar  function 

det[ln+P(s)C(s)F(s)].  tc(s)tp(s)*£(s)  (47) 

is  devoid  of  zeros  in  Re  s  >  0. 

Proof.  Theorem  2  and  (46),  Q.  E.  D. 

In  general,  independently  of  any  considerations  of  controllability  and  observ¬ 
ability, 

Ac(s)  =hc(sHc(s)  , 

Ap(s)  =hp(s)tp(s)  ,  (47a) 

Af  (s)  =  hf(s)i)f(s)  , 


Physically  speaking,  this  means  that  the  stability  of  a  syBtem  (an  internal  notion) 
cannot  always  be  predicted  from  an  examination  of  its  several  describing  transfer 
matrices  (an  external  notion). 

The  characteristic  denominator  of  a  rational  matrix  A(s)  (square  or  otherwise)  is  by 
definition  the  monic  least  common  multiple  of  all  denominators  of  aU  minors  of  A(s). 
It  is  understood  of  course  that  each  minor  is  expressed  ae  the  ratio  of  two  relatively 
prime  polynomials. 
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where  h^s),  hp(s)  and  h^(s)  are  three  polynomials  whose  zeros  determine  the  exponential 
growth  rates  of  those  "hidden"  modes  not  detected  in  the  poles  of  C(s),  P(s)  and  F(s). 

The  correctness  of  this  assertion  is  certainly  intuitively  obvious  but  a  strict  demonstra- 
ion  is  available.  Suppose,  to  be  specific,  that  the  plant  is  not  completely  controllable 
and  observable.  Then,  there  exists  a  nonsingular  constant  square  matrix  K,  square 
matrices  F^,  F ^ •  ^33  an<*  constant  matrices  F,  y  F£j,  F£j  such  that  [13,14] 


(47b) 


where  n.+n,4n,  =v  and 
1  2  3  p 

det (si  -F  )  =  *  (s)  .  (47c) 

nl  P  P 

Thus,  from  (47b)  and  (47c), 


Ap(s)  = 


= det (si 


-F  )  =det(sl 
vp  P  n3 


■F33)d',,*1»2-F22'-drt(sl„1-Fp'=Vf,)V) 


where 


h  (s)  =  det  (si  -F ’  ).det(3l  -F  )  ,  Q.  E.  D. 
p  j  j  ii^ 

Corollary  1.  If  the  "hidden"  modes  of  controller,  plant  and  feedback  sensor,  are 
all  asymptotically  stable  then,  under  assumptions  Aj,  A^,  GFCA  is  asymptotically  stable 
if  and  only  if  the  scalar  function  (47)  is  devoid  of  zeros  in  Res>0. 

Proof.  If  the  hidden  modes  of  controller,  plant  and  feedback  sensor  are  all 

jjC 

asymptotically  stable,  the  three  polynomials  hc(s),  hp(s),  hf(s)  are  strict  Hurwitz .  Now 
use  theorem  2  and  (47a),  Q.  E.  D. 

We  are  now  in  a  position  to  describe  a  truly  significant  multivariable  extension 
of  the  classical  Nyquist  criterion  for  single-input,  single-output  feedback  systems. 

Theorem  4.  Let  the  controller,  plant  and  feedback  sensor  be  asymptotically 

stable  and  define 
_ 

A  polynomial  is  strict  Hurwitz  if  all  its  zeros  lie  in  Re  s  <  0. 
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G(s)  =  F(s)P(s)C(s)  (48) 

to  be  the  GFCA  rxr  open-loop  gain  matrix  [15].  Let  r,(s)  equal  the  sum  of  all  principal 
minors  of  G(s)  formed  with  i  rows  and  columns,  (i=  1,  2. . . . ,  r),  and  put 

r 

r(s)=yr.(s).  (49) 

i=l 

Then,  under  assumptions  A^,  A^,  GFCA  is  asymptotically  stable  if  and  only  if  the 
complex  plot  of  T* ( joo)  as  U)  varies  from  to  +“  does  not  enclose  the  point  s  =  -l  in 
the  clockwise  direction. 

Proof.  Since  controller,  plant  and  feedback  sensor  are  asymptotically  stable, 
the  three  polynomials  Ac(s),  Ap(s),  A^s)  are  strict  Hurwitz.  Thus,  from  theorem  2, 
GFCA  is  asymptotically  stable  if  and  only  if 

det  [l  +  P(s)C(s)F(s)]  =  det  [l  +  F(s)P(s)C(s)] 

is  free  of  zeros  in  He  s>0.  That  is,  if  and  only  if  det  [l  +G(s)]  has  no  zeros  in  Res>0. 

r 

However,  from  a  well  known  determina.ital  expansion, 

det[lr+G(s)]  =  l+r(s)  .  (50) 

Clearly,  because  of  assumption  A  and  the  strict  Hr-witz  character  of  A  (s),  A  (s), 

I  c  p 

A^s),  the  three  transfer  matrices  C(s),  P(s),  F(s)  are  analytic  in  Res>0,  s  =“ 
included.  Moreover, 

1  +  T  (±  j”)  =  det  ( 1  r+  Jf  Jp  Jc)  =  det  [1  r r  F(»)P(®)C(“)] *  0 

and  the  standard  Nyquist  argument  [15]is  applicable  to  P(s),  Q,  E.  D. 

Corollary  1.  Let  C(s),  P(s)  and  F(s)  be  analytic  in  Re  s  >0.  Then,  under 
assumptions  A^-A^,  GFCA  is  asymptotically  stable  if  and  only  if  the  complex  plot  of 
-(jui)  as  U)  varies  from  -*  to  +"  does  not  enclose  the  point  s  =  -1  in  the  clockwise 
direction. 

Proof.  Under  A^-A^,  the  analyticity  of  C(s),  P(s)  and  F{s)  in  Re  s>0  implies 
the  asymptotic  stability  of  controller,  plant  ar.d  feedback  sensor,  respectively.  The 
rest  now  follows  from  theorem  4,  Q.  E,  D. 

>Jc 

Of  course,  if  for  some  tl‘ ,  -cc<lu<00,  T  ( j  uu)  =  - 1 ,  GFCA  is  not  asymptotically  stable. 
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*  i  iwwnpw 


If  GFCA  is  asymptotically  stable  T(s)  is  certainly  analytic  in  Re  s  >0  (Eq.  (28a)) 
but  the  converse  is  not  necessarily  true  .  Nevertheless,  in  one  special  but  very  practical 
situation  this  converse  is  valid. 

Corollary  2.  Suppose  controller,  plant  and  feedback  sensor  are  individually 
asymptotically  stable  and  assumptions  A^,  attain.  Then  GFCA  is  asymptotically 
stable  A  and  only  if  T(s)  is  analytic  in  Re  s  >0. 

Proof.  The  "only  if"  part  is  obvious.  From  the  assumptions,  A^(s),  £p(s)  and 
Aj(s}  are  strict  Hurwitz  whence,  C(s),  P(s)  and  F(s)  are  analytic  in  Re  £^0.  Suppose 
therefore  that  T(s)  is  analytic  in  Re  s^O.  Then 


T(s)F(s)=l  -[1  +P(s)C(s)F(s)]‘1 


is  also  analytic  in  Res^O  which  implies  det  [1  +  P(s)C(s)F(s)]^  0,  Res>0,  and  GFCA 
is  asymptotically  stable  (theorem  4),  Q.  E.  D. 

Under  the  conditions  of  theorem  4  an  experimental  procedure  for  measuring 
G(juj)  is  readily  described.  Referring  to  Fig.  2  observe  the  break  in  the  system  between 
the  comparator  and  controller.  With  the  input  u(t)  shorted  apply  to  the  r  terminals 


Fig.  2.  Measurement  scheme  for  G(juu),  the  open-loop  gain  matrix. 

1*,  2;, .  .  . ,  r/  the  r  steady  state  drivers  e^,  0,  0,  respectively.  Measure  the  r 

*>!< 

steady  state  responses  at  the  corresponding  points  1  ,  2  ,  .  . .  ,  r  and  amalgamate  them 
into  a  column-vector  a^(juu)e^t.  Setting 

—  - 
A  counterexample  follows  shortly. 

Because  of  the  assumptions  these  steady-state  responses  exist. 
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f 


Xi<i*>  =  •  * 


we  have 


ajOa))  =-G(jo))v1(ju)}  . 

Similarly,  with  terminal  k'excited  by  and  the  others  shorted,  let  &.  (jwje^*1*  denote 
the  measured  column-vector  steady- state  response  at  points  1  ,2  , . . .  ,  r  and  put 


vk(ju.)  =  1  I  ,  (k=  1,  2, .  , . ,  r) 


Obviously  if  R(ju;)  =  a^jui). ....  ar(jU))]» 

G(jU))  = -R{jai)  (51) 

and  the  measurement  of  G(ju))  has  been  accomplished.  In  the  practical  applications 
C(s).  P(s),  F(s)  and  therefore  T'(s)  are  functions  of  one  or  several  parameters  and  it 
becomes  necessary  to  study  the  complex  plot  of  r(jui)  over  the  range  of  these  parameters 
and  it  is  perhaps  not  out  of  place  here  to  mention  that  some  interesting  investi|ations 
along  these  lines  have  already  been  carried  out  by  Curtiss  [15]. 

A  specialization  of  theorem  3  to  the  single-input,  single-output  case  may  not 
be  without  interest.  For  this  situation  r  =n=  1,  F(s)  is  lxl  (a  scalar),  C(s)  is  mxl, 
P(s)  is  1  xm  and  P(s)C(s)  is  1  xl  (a  scalar).  Let 

a>> 

C|'|snir'  <52) 

c 


*  th 

v^  has  a  1  in  the  k —  row  and  zeros  everywhere  else. 

Transmission  from  the  primed  to  the  unprimed  side  of  the  break  is  governed  by  the 
transfer  matrix  -F(s)P(s)C(s)  =  R(s)  which  is  sometimes  called  the  return-ratio  matrix. 


14 


(53) 


and 


a>) 

p(s)  =  i^r’ 


0(s)  =ap{s)ac(s) 


af(s) 
F(fl)  =  T"7TT 


KfV* 


(54) 


where  ac{s),  a^s)  are  two  m*l  pohmomial  vectors  and  gc(s),  g^s),  gj(s),  a^s),  0(s), 
are  five  scalar  polynomials.  Subject  to  the  understanding  that  corresponding  numerators 
and  denominators  are  relatively  prime,  assumptions  A^-A^  are  satisfied  if  and  only  if 

1.  no  entry  in  a^s)  has  degree  exceeding  degree  gc(s),  n°  entry  in 
a^s)  has  degree  exceeding  g^(s)  and  no  entry  in  a^(s)  has  degree 
exceeding  degree  g^(s). 


2.  limit 


and 


0(s)  af(s) 
gjs)gjs)g.(s) 


t  -1 


3-  *c.p,£<s|'-Sc,p,f<s)- 

Theorem  3,  Corollary  2.  A  single-input,  single-output  GFCA  satisfying  assump¬ 
tions  A^-A^  is  asymptotically  stable  if  and  only  if  the  polynomial 


gp(®)gc(s)gf(s)  +  fi(®)af(s) 


(55) 


is  strict  Hurwitz. 

Proof.  Theorem  3  and  the  equalities 

Sc,p,f(s)  =Ac,p,f(s)’  Q-E-°‘ 
Observe  that  the  zeros  of  (55)  do  not  coincide  with  those  of 

0(s)a,(s) 


1  ^ 


gc(s)gp(s)gf(s) 


(56) 


unless  0(s)aj(s)  is  relatively  prime  to  g;c( s) g  ( ic )g^( s).  Thus  stability  cannot  always  be 
checked  by  simply  locating  the  zeros  of  1  +  P(s)C(s)F(s)  or  the  poles  of 
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T(s) 


P(g)C( s) 

s)P(S) 


As  an  example  suppose  c?(s)  is  a  strict  Hurwitz  polynomial  and  consider  the 
single-input,  single-output  arrangement  of  Fig.  3.  To  apply  (5?.)  we  must  examine  the 


Fig.  3.  Externally  stable  but  internally  unstable  feedback  arrangement, 
roots  of  the  equation 


or 


cp(sh)(-8)  +kiv>{-s)  =  0  ; 
cp(-s)  =  0  ;  cp(  s)  +  k  =  0  . 


Obviously,  all  the  zeros  of  co(-s)  are  in  Re  s  >0  and  the  feedback  circuit  of  Fig.  2  has 
modes  which  blow  up  exponentially  as  t-*  00  irrespective  of  the  gain  setting  k.  However 


1  +  P(b)F(«)=SI*H± 

and 

Hence,  for  sufficiently  small  |k( ,  1  +P(s)F{s)  has  no  zeros  and  T(s)  no  poles  in  Re  s>0. 

To  illustrate  the  use  of  theorem  3  we  choose  r  =  n  =  2,  anity  controller  (C  =  1^), 

,  k  real, 


\ .  s-1 

1  ' 

P(s)  = 

k  s-2 

s+2 

2 

16 


and 


F(8)  = 


fl+2s 

3s+l  ' 

i+s 

2(l+s) 

L  1 

1 

Obviously  , 


tp(s)  =  (s-2)(s+2)  -  s  -4  , 


tf(s)  =  s  +  l 


and 


det[l2+P(«)F(®)]=det 


2k+l 

~5~ 


=  |(k  +  2) 


/.  ri  =  det[l2+P(»)F(»)]^  0  ifk^-2 


By  straightforward  calculation, 

r|Ao(8)=det  [12+  P(s)F(8)]t(s}^(s)=(s+2)(3s+2)[(2+k)s-(4+k)]. 


The  zeros  of  A  (s)  are 
o 


The  latter  lies  in  Re  s<0  only  if 


p  ¥f 


i  iii  4+k 
s--2,  -2/3, 


-4  <k  <  -2 


(57) 


which  is  therefore  the  stability  range  for  the  parameter  k  . 

3.  Controllability  and  Observability  of  GFCA. 

As  we  have  already  seen  the  asymptotic  stability  of  GFCA  does  not  always  follow 
from  the  analyticity  of  T(s)  in  Re  s>0  even  under  assumption  A^.  Generally  speaking, 
the  reason  for  this  failure  is  due  to  the  fact  that  complete  controllability  and  observability 
of  controller,  plant  and  feedback  sensor  do  not  always  guarantee  that  of  GFCA  when 
viewed  as  a  system  with  input  u(t)  and  output  y(t).  Thus  it  is  possible  for  GFCA  to 


det  F(s)  h  1/2. 

i* 

It  has  been  assumed  in  accordance  with  the  requirements  of  theorem  3,  that  plan*  ^nd 
feedback  sensor  are  completely  controllable  and  observable. 
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ft 

ti 


possess  unstable  inodes  of  x  (t)  which  are  either  uncontrollable,  unobservable  or  both 
and  are  not  detectable  in  T(s).  To  examine  this  matter  in  depth  let  us  rephrase 
assumption  A^  in  terms  of  McMillan  degrees  . 

Assumption  3  (A^). 

vc=MC)  , 

Vp  =  6(P)  ,  (58) 

vf  =  6(F)  . 

In  view  of  this  assumption,  GFCA  constitutes  a  realization  of  T(s)  employing  exactly 
6(C)  +  6(P)  +  6(F)  energy  storage  elements  and  since  a  minimum  of  6(T)  energy  storage 
elements  is  required  in  any  realization  of  T(s)  [4,  5,  6,  11]  it  must  follow  that 


6(T) <  6(C)  +  6(F)  +  6(F)  . 
A  direct  proof  of  (59)  is  easily  inferred  from  the  identity 


’1 

n 

PC  ' 

—  A 

( 1  +  PCF)”  * 

-T 

* 

1 

r 

-F(  1  +PCF)" 1 
L  n 

(1  +FPC)"1 

(59) 


(60) 


Consequently  , 


taking  degrees  of  both  sides, 


6(1:<  5(PC)  +  6(F)  <  6(C)  +  6(P)  +  6(F),  Q.  E.  D. 


T(s)  is  determined  solely  by  the  completely  controllable  and  observable  part  of 
GFCA  [5,  6,  7]. 

The  McMillan  degree  of  a  ratio  ial  matrix  A(s)  is  denoted  by  6(A)  <  £  6[A(s)].  The 
reader  is  assumed  to  have  a  good  understanding  of  this  important  concept  and  a  study 
of  one  of  the  accessible  treatments  in  Refs.  4,  8- 11  is  strongly  advised. 

96c  _  j 

For  any  bv’O  rational  matrices  A(s)  and  B(s),  6(A)  =  6(A  ),  6(AB)<  6(A)  +  6(B), 

6(A  +  B)  =  6(A)  +  6(B)  and  6(K)  =  0  if  and  only  if  K  is  a  constant  matrix.  Furthermore, 
the  degree  of  a  matrix  is  greater  or  equal  to  the  degree  of  any  one  of  its  submatrices. 
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It  is  not  difficult  to  ahow  that  subject  to  assumptions  A^-A^  the  condition 

6(FPC)  =  6(PCF)  =  6(C)  +  6(P)  +  6(F)  (61) 


is  sufficient  for  the  complete  controllability  and  observability  of  GFCA  when  viewed  as 

$ 

a  system  with  input  u(t)  and  output  y(t). 

Proof.  Considered  as  a  system  with  input  u(t)  and  output  e(t),  GFCA  is  described 
by  the  transfer  matrix  (16): 

E(s)  =  [lr+F(s)P{s  JQs)]"1  . 


Since  6(E)  =  6(FPC), 


6(FPC)  =  6(C)  +  6(P)  +  6(F)  (62) 

certainly  suffices  to  guarantee  the  controllability  of  GFCA  with  input  u(t)  and  output 
y(t).  By  duality  [6],  the  observability  of  GFCA  with  input  u(t)  and  output  y(t)  is  equivalent 
to  the  controllability  of  the  dual  dynamical  system 


Sd(t)*Fo2d(t)  +  HoHd(t)  * 

(63) 

id(t)=Goid(t)  +  JoHd(t)  • 

(64) 

The  reader  should  have  no  trouble  in  convincing  himself  that  the  dual  GFCA,  (GFCA)^ 
say,  is  correctly  represented  in  Fig.  4  provided 

C^-fF'.H^G'.j;), 

P'(s)-(F',H'  G'  j')  , 

P  P  P  P 

F/(s)-(Ff,,H',Gj,jJ) 

where 

C(s)  **(F  ,  G  ,  H  ,  J  )  , 
c  c  c  c 

P(s)  -*(F  ,  G  ,  H  ,  J  )  , 

P  P  P  P 

F(p)-(F{>  Gf,Hf,  J£)  . 

* 

This  is  lemma  1  of  Chen's  paper  [2].  In  a  private  communication  to  the  author  Chen 
has  outlined  a  different  proof  based  on  some  earlier  results  obtained  in  Ref.  1.  Clearly, 
as  far  as  controllability  is  concerned  the  choice  of  output  point  is  immaterial. 
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Thus  C'(a),  P'(s),  F#(s)  are  all  realised  minimally  and  is  therefore  also  satisfied  by 
(GFCA)d. 


Fig.  4.  A  schematic  cf  (GFCA).,  the  generalized  feedl  '<  control  arrangement  dual 
to  that  of  Fig.  1. 

Viewed  as  a  system  with  input  u^tt)  and  output  e^t),  (GFCA)^  has  the  transfer 
matrix  description 

E,(s)  =[1  +F/(s)C/(s)P,(s)]‘1.  (65) 

d  n 

Hence  controllability  of  (GFCA)^  and  concomitantly,  the  observability  of  GFCA  with 
output  y(t)  is  assured  if 

SCF'C'P')  =  6(C')  +  6(P')  +  6(F'); 
or,  since  a  matrix  and  its  transpose  have  the  same  degree,  if 

6(PCF)  =  6(C)  +  6{P)  +  6(F)  .  (66) 

Taking  (62)  and  (66)  together  we  get  (61),  Q.  E.  D. 

Instead  of  working  with  the  criterion  (6 1 )  which  is  rather  intractable  anyway 
it  appears  advantageous  at  this  stage  to  introduce  four  realistic  and  simplifying  restric¬ 
tions  from  the  very  outset. 

Rj.  The  number  of  plant  inputs  m  equals  or  exceeds  the  number  of  plant 
outputs  n  ;  i.  e. ,  m  >  n. 

That  R.  is  physically  correct  is  obvious  and  that  it  is  satisfied  in  almost  all  practical 

*  sjc 

situations  is  attested  to  by  the  abundance  of  examples  described  in  the  literature  . 

R^.  The  number  of  system  input  variables  equals  the  number  of  output 
variables  n;  i.  e.  ,  u(t)  and  y(t)  have  the  same  dimension  and  r  =n. 

See,  for  example,  the  interesting  book  by  Meerov[l2], 
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Since  u(t)  serves  as  the  reference  vector  or  "cetpoint"  for  the  plant  output  y(t),  the 
plausibility  of  should  be  apparent. 

Rj.  The  plant  transfer  matrix  P(s)  is  non- degene  rate.  That  is,  P(s) 
has  normal  rank  equal  to  the  number  of  its  rows  : 

normal  rank  P  =  n  . 

If  R^  is  false  there  exists  at  least  one  polynomial  dependence  between  the  rows  of  P(s) 
This  implies  in  turn  the  existence  of  a  fixed  nontrivial  set  of  polynomials 
a^(s),  a^s), . . .  ,  an(s),  independent  of  u(s),  such  that 


a1(s)y1(s)  +  a2(s)y2(s)  +  . . .  +  an(s)yn(s)  s  0  ,  (67) 

where  the  jMs)  are  the  components  of  y(s),  (i  =  1,  2, . .  . ,  n).  Hence,  at  least  one  trans¬ 
formed  output  variable  may  be  eliminated  in  terms  of  the  others  and  therefore  becomes 
superfluous. 

R^.  The  mxn  controller  transfer  matrix  C(s)  is  chosen  so  that  the  square 
nxn  matrix  P(t-)C(s)  has  normal  rank  n;  i.  e. , 

det  [P(s)C(a)]  ^  0  .  (68) 

Choosing  C(s)  so  that  det  [P(s)C(s)]  aO  leads  to  an  identically  singular  T(s): 

detT(s)  *  0  .  (69) 

Again  this  implies  a  polynomial  dependence  of  the  type  (67)  on  the  n  transformed  output 

variables  and  GFCA  is  deficient  from  the  point  of  view  of  output  controllability.  Stated 

differently,  under  (69)  there  exists  a  very  rich  class  of  reference  inputs  (including  even 

some  of  the  "step"  type)  which  GFCA  cannot  follow.  These  practical  considerations 

coupled  with  the  relatively  simple  theory  that  results,  is,  in  our  opinion,  ample 

justification  for  pursuing  the  consequences  of  a  design  procedure  modeled  around 

R,  -R .. 

1  4 

Definition  1.  A  GFCA  satisfying  assumptions  A^-A^  and  itrictions  R^-R^ 
is  said  to  be  standard. 

Theorem  5.  A  standard  GFCA  is  completely  controllable  and  observable  if  and 

only  if 

& 

A  rational  matrix  A(s)  is  said  to  have  normal  rank  equal  to  k  if  1)  at  least  one  kxk 
minor  does  not  vanish  identically  but  2)  all  minors  of  orders  larger  than  k  do  vanish 
identically. 
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1.  6[(PC)’ l  + F]=  6[{PC)’1  ]  +  6(F)  and 

2.  6(PC)  =  6(P)  +  6(C). 

Proof.  Invoking  (68), 

T=PC(1  +FPC)"1=*[(PC)"1+F]~1  . 

.’.  6(T)  =  6[(PC)-1+  F]<6[(PC)-1]+  6(F)  =  6(PC)  +  6(F)  <  6(C)  +  6(P)  +  5(F)  . 

(70) 

Thus,  to  achieve  equality  in  (70)  between  the  extreme  left  and  right  members  it  is 
necessary  and  sufficient  to  have  both 

6[(PCf  1+F]=6[(PC)'1]+6(F) 

and 

6  (PC)  =  6(P)  +  6(C),  Q.  E.  D. 

Corollary  1.  A  sufficient  set  of  conditions  for  a  standard  GFCA  to  be  completely 
controllable  and  observable  are 

S^.  (PC)-1  and  F(s)  have  no  common  poles. 

32.  5(PC)  =  6(P)  +  5(C). 

In  the  single-input,  single-output  case  (r  =n  =  l)  and  are  also  necessary. 

Proof.  It  is  a  property  of  McMillan  degree  that  6(A+ B)  =  6(A)  +  6(B)  whenever 

A(s)  and  B(s)  are  devoid  of  common  poles,  Now  let  r=n  =  l,  Then  C(s)  is  mxl, 

P(s)  is  1  xm,  PC  is  lxl,  F(s)  is  lxl.  Thus  (PC)"^  and  F(s)  are  two  scalar  rational 

functions  a(s)  and  b(s),  say.  If  a(s)  and  b(s)  admit  s  =  sq  as  a  common  pole  of  respective 

orders  r,  >1  and  r,>l,  it  is  clear  that  s  =  s  ,  as  a  pole  of  a(s)  +  b(s),  cannot  have  order 
i  c  o 

exceeding 


max(rj,  r,)  <  r^  r2 

and  6(a+b)  =  5(a)  +  5(b)  is  impossible,  Q.  E.  D, 

Corol.  ary  2.  A  standard  GFCA  with  unity  controller  (C(s)  =  1^;  n  =  m  =  r)  is 
completely  controllable  and  obs  ;rvnble  if  and  only  if 

6(P_1+F)  =  6(P_1)  +  5(F)  .  (71) 

— 

And  so  we  have  the  classical  rule  forbidding  the  cancellation  of  any  zero  of  the  net 
feedforward  transfer  function  by  a  pole  of  the  net  feedbackward  transfer  func^'on. 
Incidentally,  this  shows  that  (61)  is  not  even  necessary  in  the  single-input,  single¬ 
output  case! 
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Moreover,  to  achieve  equality  in  (72)  it  suffices  that  P~*(s)  and  F(s)  be  devoid  of 
common  poles. 

Despite  the  extensive  analysis  presented  above  the  constraint  of  complete 
controllability  and  observability  is  not  necessarily  an  important  nor  even  desirable 
attribute  of  a  GFCA!  This  statement  is  undoubtedly  disconcerting  and  puzzling  but 
consider  that  the  main  objective  of  GFCA  is  not  to  control  or  observe  the  internal 
state  x  (t)  but  to  drive  the  plant  output  y(t)  to  the  "set  point"  u(t).  And  in  fact  it  seems 
safe  to  say  that  most  successful  multivariable  control  systems  are  not  designed  under 
this  constraint  which,  when  superimposed  on  the  other  requirements  of  steady- state 
performance,  transient  response  and  asymptotic  stability,  presents  insurmountable 
algebraic  difficulties.  What  is  of  paramount  importance  is  that  GFCA  be  asymptotically 
stable.  As  we  know  from  theorem  4,  corollary  2,  there  are  many  practical  situations 
in  which  this  is  ascertainable  purely  from  an  examination  of  the  analyticity  properties 
of  T(s)  in  ResiO  even  if  GFCA  is  not  completely  controllable  and  observable.  In  any 
case  theorems  2-4  and  their  accompanying  corollaries  provide  an  adequate  apparatus 
for  checking  the  stability  of  GFCA  irrespective  of  other  considerations.  Also,  as  is 
well  known,  the  use  of  redundant  energy  storage  elements  may  improve  both  the 
sensitivity  and  reliability  of  GFCA,  permit  the  use  of  more  realistic  topologies  and  even 
facilitate  the  synthesis  problem  of  minimizing  one  or  several  cost  indices.  Nevertheless 
the  effort  spent  in  this  section  is  not  wasted  because  it  has  revealed  insights  which  not 
only  deepen  our  understanding  but  also  provide  the  conceptual  background  essential  to 
the  development  of  perspective. 

5.  The  Problem  of  Plant  Equalization.  The  given  datum  is  the  nxm  transfer 

* 

matrix  P(s)  of  a  dynamical  plant  with  asymptotically  stable  hidden  modes  . 

Definition  2,  The  nxr  matrix  T(s)  is  said  to  be  realizable  for  P{s)  if  for  some 
choice  of  asymptotically  stable  dynamical  controller  and  feedback  sensor  the  resultant 
GFCA  of  Fig.  1  incorporating  the  given  plant  is  a  dynamical  asymptotically  stable  system 
possessing  the  prescribed  input-output  transfer  description  T(s). 

Analytically  the  question  of  realizability  reduces  to'the  following:  Given  a 
P(s)  which  is  finite  at  s=“  and  a  T(s),  find,  if  possible,  two  real  rational  matrices 


Thus,  the  stability  properties  of  the  plant  are  completely  summarized  by  ^(s),  the 
characteristic  denominator  of  P(s)  (See  theorem  3,  corollary  1). 
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$ 

C(s)  and  F(s)  both  analytic  in  Re  s>0,  s=®  included,  such  that 

1.  T(s)  =  Cln+P(s)C{s)F(8)]-1P(s)C(s)  ; 

2.  det  [1r+  P{®)C(®)F{»)]  t  0  and 

3.  det[ln+P(s)C(s)F(s)]ys)  ^  0,  Re  s>0  .  (72) 


It  is  intended  that  C(s)  and  F(s)  shall  be  realized  minimally;  i.  e.  ,  as  completely 
controllable  and  observable  dynamical  systems.  Thus  A  (s)=s}[  (s)  and  a  ( s }  =  ih( s ) 
are  strict  Hurwitz  and  3  reduces  to  the  stability  criterion  of  theorem  2.  Part  II 
(in  preparation)  presents  a  complete  solution  for  the  standard  GFCA  together  with  a 
full  discussion  of  its  theoretical  and  practical  implications. 
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